Two-Point Versus Multipartite Entanglement in Quantum Phase Transitions 
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We analyze correlations between subsystems for an extended Hubbard model exactly solvable in one dimen- 
sion, which exhibits a rich structure of quantum phase transitions (QPTs). The T = phase diagram is exactly 
reproduced by studying singularities of single-site entanglement. It is shown how comparison of the latter quan- 
tity and quantum mutual information allows one to recognize whether two-point or shared quantum correlations 
are responsible for each of the occurring QPTs. The method works in principle for any number D of degrees of 
freedom per site. As a by-product, we are providing a benchmark for direct measures of bipartite entanglement; 
in particular, here we discuss the role of negativity at the transition. 
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In the past years, the characterization of complex quan- 
tum phenomena has received a strong impulse from the recent 
developments in quantum-information theory. Within such 
framework, a crucial notion is that of entanglement. Besides 
being recognized as a fundamental resource for quantum com- 
putation and communication tasks 1 1], it has also been used 
to better characterize the critical behavior of different many- 
body quantum systems when some characteristic parameter of 
the related Hamiltonian is varied; the latter phenomenon be- 
ing known as quantum phase transition (QPT) 1 2]. 
In fact, a deep comprehension of universal properties of QPTs 
has not been fully reached yet. The peculiarity of using entan- 
glement in this context is that, being a single direct measure of 
quantum correlations, it should allow for a unified treatment 
of QPTs; at least, whenever the occurring QPT is to ascribe to 
the quantum nature of the system, which is always the case at 
T = since thermal fluctuations are absent. 
A first description of the relations between entanglement of 
one or two spins and QPTs in spin- 1/2 chains was given in 
01, where it was noticed how derivatives of concurrence show 
divergencies in correspondence of QPT, with appropriate scal- 
ing exponents. The entanglement of blocks of L spins and its 
scaling behavior in spin models showing critical behavior was 
then investigated in |4|]. The problem of characterizing the 
ground state phase diagram of fermionic systems by means of 
entanglement has been addressed more recently in |5], where 
it was shown how the study of single-site entanglement allows 
one to reproduce the relevant features of the known (numeri- 
cal) phase diagram. While this is a promising starting point, it 
remains to be clarified which quantum correlations are respon- 
sible for the occurring QPT: if two-points or shared (multipar- 
tite), if short or long ranged. The answer to the above issue 
would in fact require exhaustive investigation of the entangle- 
ment between any two subsystems. In case the subsystems 
have just two degrees of freedom, concurrence properly quan- 
tifies the quantum correlations Iffl. A generalization of such 
quantity to (sub)systems with a higher number of degrees of 
freedom D has been proposed, and is known as negativity |7]. 
Also, the total amount of correlations between any two sub- 
systems is captured by quantum mutual information 0]. 



In the following we describe a method based on the compar- 
ison of the latter quantities for arguing whether the occurring 
transition is to ascribe to two-points or multipartite quantum 
correlations; the method works for arbitrary D. Our strategy 
is tested on a one-dimensional extended Hubbard model that 
was solved |3, llfJl . exhibiting a rich structure of phase dia- 
gram at T = 0. We show that the phase diagram is exactly re- 
produced by the singularities of single-site entanglement. We 
then infer which of the QPTs is originated from a singular be- 
havior of two-point or multipartite entanglement; our results 
are confirmed by the exact solution. 

Correlations and subsystems. We are interested in the ex- 
isting correlations between: a) the single site i and the rest of 
the system; b) the generic site i and a generic site j ^ i;c) the 
generic pair of site (dimer) and the rest of the system. 
In order to measure the total correlations between two generic 
subsystems A and B we use the quantum mutual information 
1 lL IsL 1 1 ill . The latter is defined as 

Tab = S{p A ) + S{p B ) - S(pab), (1) 

where pab, Pa and ps are the system's and subsystems' den- 
sity matrices, respectively, and S(p) = — J^i \ l°g2 (^i 
being the eigenvalues of p) is the Von Neumann entropy. In 
O llltl it was shown how Tab is a proper measure of all 
(quantum and classical) correlations between A and B. In 
case A and B are single sites, we will refer to the latter as 
two-point quantum (Q2) and classical (C2) correlations. 
As far as quantum correlations are concerned, we consider 
two different cases. When pab is a pure state, correlations 
between A and B are purely quantum and are measured by 
S(pa) = S(pb)- This happens when A corresponds to a sin- 
gle site % (or to the dimer and B corresponds to the re- 
maining sites 11211 : Si = S(pi) (single-site entanglement) ac- 
counts for both the localized correlations (Q2) and the shared 
ones (QS in the following). When we deal instead with the 
correlations between two generic sites the density ma- 

trix of the global system is the dimer's one: pab — Pij- The 
latter generally corresponds to a mixed state. Thus, to evalu- 
ate the quantum correlations between two generic sites, we 
need a measure of entanglement for bipartite mixed states. 



In general, proposed measures are hard to compute whenever 
D > 2, since they require difficult optimization processes. 
However, there is at least one measure easy to compute 0, 
the negativity 



^(pa b ) = (||p^||i-1)/2; 



(2) 



where p/^ is the partial transposition with respect to the 
subsystem A applied on pab, and ||0||i = TrvOH) is 
the trace norm of the operator O. p T ^ B can have negative 
eigenvalues pi, and the negativity can also be expressed as 
N(pab) — | Yli Mil- Although negativity is not a perfect 
measure of entanglement 1 13], it gives important bounds for 
quantum information protocols i.e., teleportation capacity and 
asymptotic distillability. Its role in describing QPTs has not 
been fully investigated yet. 

Entanglement and QPTs. Si has been proven to be a use- 
ful tool in describing QPTs |5]. As already pointed out, to 
give a better characterization of the latter, one could as well 
consider quantum correlations between different subsystems. 
The scheme we propose in this letter is based on the idea of 
comparing Si -not allowing one to distinguish Q2 from QS 
correlations- with different functionals quantifying instead 
just two-point correlations. We study Afi.j, which is at least 
a lower bound for Q2 correlations, and li.j, which properly 
captures total (Q2 and C2) correlations. As a first step, the 
exact phase diagram is obtained analyzing the singularities 
shown by Si, Tij, and Afij. Successively, a comparison of 
the singular behavior of Zj j with that of Si allows one to dis- 
criminate whether a QPT is to ascribe to Q2 or QS correla- 
tions. In fact, whenever Si exhibits a singular behavior due to 
Q2 correlations, the same type of singular behavior should be 
highlighted as well by Ti j (since it also contains Q2 correla- 
tions), and possibly by Afi.j, in case the latter would properly 
capture them for our model. On the contrary, when the singu- 
lar behavior of Si is to ascribe to QS correlations, the same 
singular behavior should not be displayed either by Zjj or by 
Ni,j, since both measures regard only two-point correlations. 

The bond-charge extended Hubbard model. The model 
we deal with is described by the following Hamiltonian: 



H 



BC 



[l-xinw+nj^jcl^Cja , (3) 



<i ,j><7 



where c i(T , c ia are fermionic creation and annihilation opera- 
tors on a one-dimensional chain of length L; a =|, | is the 
spin label, a denotes its opposite, n JCT = c^c^ is the spin 
a electron charge, and (i, j) stands for neighboring sites; u 
and x (0 < x < 1) are the (dimensionless) on-site Coulomb 
repulsion and bond-charge interaction parameters. 
The model is considered here at x = 1, in which case the 
number of doubly occupied sites becomes a conserved quan- 
tity. The exact eigenstates of Hbc{x — 1) are obtained in 
||9j, [lOj Il4ll : the ground-state phase diagram is shown in the 
left part of Fig. \l\ The latter presents various QPTs driven by 
parameters u and average number of electrons per site (filling) 
n. The charge-gapped phase IV is insulating and all sites are 
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FIG. 1: Left: Ground-state phase diagram. Empty, slashed and full 
dots stay for empty, singly and doubly occupied sites. Right: d u Si. 



singly occupied; phases I, and 77 fall in the Tomonaga- 
Luttinger class (neither spin nor charge gap); they are charac- 
terized by the presence of singly and empty sites (phase I), 
singly and doubly occupied sites (phase I'), both of which 
have dominant charge-charge correlations, and all types of 
sites (phase II) with superconducting correlations and off- 
diagonal long range order (ODLRO). The latter characterizes 
also phase 177, where sites are empty or doubly occupied. 
The model's energy spectrum is fully independent of spin ori- 
entation 1 9]: any sequence of spins in the chain cannot be 
altered by the Hamiltonian, which, in fact, acts on a Hilbert 
space that at each site i has Di = 3, and is spanned by the 
states 1 0) / (empty), \0)i (singly occupied), and \2)i (doubly 
occupied). The physics of the system is essentially that of 
N s = Y^i cjcj spinless fermions and Nd bosons, with eigen- 
states given by 



(4) 



Here A/" 



1/2 



(L-N s - N d )\/{L - N s )\N d \ 
malization factor; a}- is the Fourier transform of the spin- 
less fermion operator cj, a] = J2 q ex P(^x^) c 9' 
j = 0, . . . , N s — 1; moreover 7/^ = 5D i=1 T]J is also known 
as the eta operator, and creates doubly occupied sites from 
empty ones (f^lO)^ = |2)i); \vac) is the electron vacuum. 
{rp ) k |?;ac]_is known to carry ODLRO and multipartite entan- 
glement 0]. At fixed filling n = (N s + 2N d )/L, the actual 
value of N s in @ is chosen to minimize the corresponding 
eigenvalue E(N S , N d ) = -2- sin(Tr^) + uN d . 



The system density matrix in the ground state is defined by 
p = \ip(N Sl N d ) >< ip{N s ,N d )\. Results of the calculation 
for the single-site pi and the dimer reduced density matri- 
ces are reported below. With respect to the basis |0), |0), |2), 
Pi = diag{l — n s — n d , n s ,n d } with n a = ^ (a = s, d). 
Whereas with respect to the basis |00), |O0), |0O), |00), |02), 
1 20), |02), |20), |22), 
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Here, assuming e = 1/L, 
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In the thermodynamic limit e — > 0, n a finite, the above results 
may also be derived from 1 10]. 

Results. As a preliminary observation let us notice that in 
phases I, I' and III (see Fig. ^ left side) the dimension 
of on-site vector space reduces to two, meaning that in these 
cases JVi.j should reproduce results evaluated through concur- 
rence. This happens to be the case; in particular, in phase 77/ 
Mi j (and the concurrence) are vanishing V|i — j\, whereas 
Ti.j is equal to n{2 — n)/2, which is related to the value of the 
ODLRO parameter, in agreement with lfl5ll . 
We also observe that whenever Cy is zero (for instance, 
phases III and IV) the two-site density matrix is indepen- 
dent of the sites i and j. In the insulating phase IV this hap- 
pens because the state is a tensor product of identical single- 
site states, and all correlations are identically zero. On the 
contrary, in phase III, since Q ^ 0, the surviving two-point 
quantum correlations are range independent. In such cases it 
may be useful to introduce global (i.e. sums over all sites) 
quantities instead of local ones, since it may happen that a 
correlation is locally vanishing but globally relevant; also fi- 
nite size corrections (e) have to be considered. For instance, 
it turns out that in so doing in phase III the total negativity 
becomes nonvanishing = "(2 — n)/[(2— n) 2 +n 2 ]. 

We now turn to discuss what happens at QPTs by study- 
ing the behavior of Si, Ti.j, Mi.j. As mentioned, each of the 
observed measures of correlations keeps track of the under- 
going transitions, exhibiting a singular behavior at the tran- 
sition points. The latter can be characterized by the analysis 
of the partial derivatives of each measure. As an example, in 
the right part of Fig. ^we plot d u Si. Noticeably, the diver- 
gencies in the derivative are in perfect correspondence with 
the parameter's values at which the various QPTs occur, aside 
from transitions I, I' — ► IV that must be revealed by d n Si. 
The systematic analysis of the behavior of the various deriva- 
tives at each QPT is carried on in Table H] We first consider 
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TABLE I: Behavior of the evaluated partial derivatives at critical 
points for the various QPTs (left column): FD is finite discontinu- 
ity, "multi" refers to multipartite, and "two" refers to two-point. 



d x Si, x — n,u; it exhibits two different kinds of divergen- 
cies: logarithmic for transitions I, I' — > IV and II — > I, I'; 
algebraic for transitions II — > 77/ and II — > IV, with expo- 
nent v = 1/2. The latter turns out to correspond to the shift 
exponent as extracted from finite size analysis 1 1(J. 
In Tab. [Qwe also report the behavior of d x Xij and d x M^j 
at QPTs. As described in the paragraph Entanglement and 
QPTs, the comparison of the three quantities can be used to 
understand whether bipartite or multipartite entanglement is 
relevant to the various transitions. In fact, all transitions cor- 
responding to a (logarithmic) divergence in d x Si are not seen 
as divergencies either in d x Iij or in d x Ni,j ■ In such cases, we 
infer that the transitions are to ascribe to QS correlations; this 
is also in agreement with the fact in some of these transitions 
(II — > 1,1') the component of the ground state given by the 
eta pairs (which carry multipartite entanglement and ODLRO) 
disappears. On the contrary, whenever the divergent behavior 
exhibited by d x Si is also displayed by <3 x 2y, as seen for the 
two transitions II — * III and II — > IV, this is to be inter- 
preted as a signal of the role of Q2 correlations in the QPT. 

Different from what happens for the first two quantities 
reported in Table |l| we can check on the third column that 
d x Nij never does display the same singular behavior as d x Si. 
In particular, this happens in correspondence of the transitions 
driven by two-point correlations, II — > III and II — > IV, 
suggesting that Afij is just a lower bound for Q2 correlations 
also for the present model. Apart from this fact, the behav- 
ior of Afi.j for various values of \i — j\ supports once more 
the idea that the transitions in question have to be ascribed to 
two-point correlations. As an example, we report in the top 
left part of Fig. HjAij for the case n = 0.5; the transition 
II — » III takes place at u c = —4. As u gets close to u c two- 
point quantum correlations begin to spread along the chain; 
this is shown by the non zero value of Mi j for an increasing 
number of pairs of sites whose distance \i — j\ grows up to 
oo as u — * u c . This is a clear indication of diverging corre- 
lation length originated from Q2 correlations at critical point. 
One could expect that again the total negativity is the right 
quantity to display a critical behavior, in agreement with sim- 
ilar conclusions about concurrence 1 17] in spin-1/2 systems. 
Moreover, the value at which Mi,j reaches its maximum gets 
closer to u c by increasing \i — j\, indicating its possible scal- 
ing behavior. The same qualitative behavior of maximum is 
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tity with singularities of single-site entanglement allows one 
to distinguish at each QPT the contribution of bipartite from 
that of multipartite entanglement. At the same time, whenever 
a contribution from two-point quantum correlations is spotted, 
this can be used to test direct measures of bipartite entangle- 
ment. As an example, we tested the negativity, finding that in 
this case it does not capture all of the two-point quantum cor- 
relations, though it shows evidence of a diverging correlation 
length and interesting scaling behavior in the vicinity of the 
transitions ascribed to two-point correlations. The study of 
scaling properties of the proposed measures of entanglement 
and of total negativity II 1611 need to be further investigated. 
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FIG. 2: Plots of A/ij, Si and Sij (\i — j\ = 1, 5) for the 
section n = 0.5 (line a in Fig.0 



observed for Zjj (top right part of figure), even though such 
a quantity is, in general, different from zero also away from 
the critical point, suggesting once more that the quantum mu- 
tual information in the vicinity of the transition captures the 
divergent behavior of just the Q2 correlations. We finally an- 
alyze in Fig. |H<Si,j = S(pij), which describes all quantum 
correlations between the dimer i, j and the rest of the system. 
Interestingly, Si.j has for all i,j the same qualitative behav- 
ior of Si at critical points. Such a feature is confirmed by our 
calculations in correspondence of all QPTs. This is expected 
within our scheme, since Si and Sij both describe the same 
correlations, Q2 and QS. 

Conclusions. We have studied the behavior of different 
measures of correlations in correspondence of QPTs for a 
solvable model of correlated electrons on a chain, displaying 
different kinds of metal-insulator-superconductor transitions. 
As a general output of our work, the role of quantum mutual 
information in the investigation of QPTs has been recognized. 
In particular, the comparison of singularities of the latter quan- 
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